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ABSTRACT 


This study investigates fin-line filter structures, a 
type of E-plane waveguide device. Expressions germane to 
the analysis of fin-line structures are developed. A 
computer aided design and analysis program based upon a 
mode-matching technique is described. Filters designed 
by this program are fabricated and tested in X-band. Good 
agreement with predicted response is obtained. Filters 
fabricated and tested in Ku and Ka bands by other researchers 
are analyzed by the new program. Good agreement between 


predicted response and published performance is noted. 
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I. INTRODUCTION 


. 


Integrated circuits for use below about 3 GHz demonstrate 
Clear advantages in terms of size, weight, cost and electri- 
cal performance, over alternative techniques. These 
advantages often accrue as a result of component miniaturi- 
Zation. However, as the frequency is increased to the 
centimeter and millimeter regimes, the resultant excessive 
miniaturization creates problems associated with critical 
mechanical tolerances, and questionable production uniformity. 
Moreover, concomitant undesireable electrical performance 
involving radiation loss, spurious coupling, dispersion and 
higher order mode propagation becomes unacceptable [Ref. 1]. 
Over the last decade, integrated fin-line structures have 
received increased attention as attractive media for 
low-insertion loss design applications which avoid excessive 
miniaturization, yet offer the potential for low-cost 
batch-processed production techniques. 

The fin-line filter structure consists simply of 
conductive strips, or fins, which bridge the broad walls of 
a rectangular waveguide. That is, for waveguide operated 
in a TE mode, the conducting fins are suspended in the 
E-plane. Metal sheet structures suspended in the center of 
the waveguide (herein referred to as "Class I" structures) 


offer the simplest possible structure for analysis and 
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design. These unilateral center-line fins present fabrica- 
tion and suspension difficulties. Bilateral fin-line filters 
may be constructed as printed circuits on low-loss dielectric 
substrates and inserted in the center of the guide. These 
("Class II") filters offer the advantage of construction by 
well developed spintéa circuit techniques, and relatively 
straight forward analysis [Refs. 2, 3], but still present 

the problem of mechanically fixing the fin-line structure in 
the center of the guide. Clearly, fin-line structures printed 
on dielectric substrate and affixed to the waveguide side- 
walls ("Class III" structures) offer mechanical advantages 
over the other two class types. The authors know of no 
previous investigations concerning Class III strucutres. 

The purpose of this thesis study was to develop and 
validate a computer aided design (CAD) program for the desig: 
and analysis of all three class types of E-plane fin-line 
filters. To this end, the appropriate mathematical 
expressions were derived for a generalized fin-line structure. 

The generalized structure has two different dielectric 
materials in the guide. The Class II filter described above 
corresponds to the case where the higher dielectric constant 
material is in the center of the guide. The Class III filter 
has the higher permittivity material along the waveguide walis 
The generalized Class I filter has the same permittivity 


throughout the fin-line structure. 
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The derived transcendental equations relating the mode 
wave number values in different zones of the waveguide 
apply to all class types. However, when the difference 
in dielectric constants is large the trigonometric functions 
in the transcendental relationship may become hyperbolic. 
Moreover, Class II and Class III structures require different 
but equivalent expressions in the numerical analysis to 
avoid singularity occurances. 

The analysis method employed is a mode-matching technique 
uniquely applied to include the effects due to an arbitrary 
number of higher order modes. A convergence dependence 
analysis assures the most efficient regional distribution of 
the limited number of higher order modes retained for the 
numerical analysis. 

Validation of the complete computer program for Class | 
I and II filters was initially conducted by analyzing Ky and 
Ko band fin-line filters designed and constructed by Shih 
feet. .-2), and Arndt, et al [Ref..31.. Finally, several Class 
I and Class III filters were designed, constructed and 
tested locally. These test filters were resonant in 


X=band (8-12 GHz) to match available test hardware. 
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Ti, THEORY 


For the types of fin-line structures in waveguides 
considered here, each geometrical discontinuity may be 
considered in two distinct coordinate frames. It is useful, 
therefore, to expand the unknown fields in both regions in 
terms of their associated normal modes. Field continuity 
requirements are then imposed on the modes at the interface 
of the regions. Finally, the orthogonality property of the 
normal modes is used to derive an infinite set of linear 
Simultaneous equations for the unknown modal coefficients. 
This is the essence of the mode-matching technique for the 
solution of boundary value problems [Ref. 4]. Moreover, 
this procedure leads directly to the scattering matrix for 
the discontinuity. Since the fin-line is a series of such 
discontinuities separated by transmission line sections, a 
cascading technique may be employed to develop an equivalent 
scattering matrix for the entire structure. 

In principle, the effects of all the possible modes upon 
each of the others is inherently included, and the analysis 
is exact. In practice, the number of modes considered must 
be limited to a reasonably finite set. Accurate and 


ient convergence of the subsequent numerical calculations 


JH 
Q 
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is achieved by following the guidelines suggested by the 


+ 


studies of Shih and Gray [Ref. 5]. 
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The derivation of the equivalent scattering matrix for 
the generalized fin-line structure follows. Detailed 
intermediate steps of the development are contained in the 
appropriate appendices. The basic geometry of the E-plane 
fin-line is shown in Fig. II-l. It is important to note 
that the assumed "magnetic wall" (x=0), along which the 
parallel magnetic field is zero, implies TE oO type (m odd) 


incident waveforms. The computer simulation is based upon a 


TEL 5 incident wave. 


A. EXPRESSIONS FOR NORMAL MODES AND EIGENVALUES 

By assuming incident modal electric (E) fields of the 
TE Lo type, field symmetries allow the simplified waveguide 
geometry shown in Fig. II-2. The conducting fin is assumed 


to be infinitely thin. The dielectric materials in the guide 
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are assumed to be lossless, homogeneous, and istropic. Al 
the fin-line structures and substrates must be constructed 
to maintain symmetry about the x = 0 plane. For simplicity 
of discussion, the guide is considered to be comprised o 
two regions, and two zones. The regions are Z dependent. 
Region I refers to the guide before the metal fin boundary, 
and Region II refers to the guide after (plus Z direction) 
the fin boundary. The zones are X dependent. Zone 1, 
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sometimes called the lower zone, refers to the condition: 
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while Zone 2, sometimes called the upper zone, refers to 


the condition: 
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Modal eigenvalues in -Region I use lower-case subscripts, 
while those of Region II use upper-case. Regional divisions 
are occassionally indicated explicitly by the superscripts. 
ee and "TT". 


The Z dependence of all fields is of the form: 


: ; acts k th 
where Yn <8 the propagation constant of the m mode. 


Likewise, the time dependence of all fields is of the form: 


where » is the radian frequency of the electromagnetic (EM) 
wave. These terms are understood to multiply all field 
expressions, and will not be explicitly shown, except where 
it is more convenient, or constructive to do so. 

The general solution to the wave equation has the form: 
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where the - are the associated wave numbers and the xX are 
the appropriate directional variables, x or y. 

For our purposes, the electric field in the Y direction, 
By is of greatest interest. In Region I, the appropriate 


expressions for the normal modes of the = field are found 


to be (Appendix A): 
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The wave numbers for the Nae mode are related by the 


transcendental equations: 
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Equation II.A.2 is the usual dispersion relation, 
where aay and e« are the relative dielectric constants of the 
two zones. The transcendental equations (II.A.2) and (II.A.3) 
are solved simultaneously to determine the value of the wave 
number for each mode in both zones. 

The corresponding expressions for Region II are similarly 


found to be: 
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The wave number for the a mode is: 
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It should be noted that in Region II, the eigenvalues for 


+ + 1 
Kip are uncoupled from the eigenvalues for Kup 


mode number, p, runs from 1 to » individually in each Zone. 


SO that the 
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The propagation number for the m~ mode in Region I is 


given by: 


Positive real values for the propagation number corresponds 
to evanescent waves. Propagating wave modes occur when 


the propagation number is imaginary. 
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B. MODAL ANALYSIS 
Any physically realizable electric field may be written 


as some linear combination of the normal modes. That is: 


Boe) x Ao. 
y n=1 
and 
= — ee Smee | ee : 
y pei 2 D 


The So and are the normal modes of the =. field in Region 
I and Region II, respectively. The A, and D_ are the 
associated amplitudes. 

The interface discontinuity may be analyzed as a two 


mart junction. with values as defined in Fig. It=}3. The 
p= a 5 So 


boundary conditions on the 9 and w_ [Ref. 4] require that: 
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G2 mths ; 23 
ain = Coefficients of Modes Incident from the Left 
A = Coefficients of Modes leaving junction to the Left 
D, = Coefficients of Modes Incident from the Right 
D = Coefficients of Modes leaving junction to the Right. 
i 
[S$] = Seattering Matrix 


Fig. II-3. Two-port Junction Model of Discontinuity 
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The characteristic impedance of the n mode, ser is defined 
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The orthogonality of the normal modes implies that: 
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where a is the kronecker delta. 


We define the coupling coefficients, H 
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By representing the A and D amplitudes as vectors, and the 
H coupling coefficients and Y impedance terms as matrices, 


these relationships may be expressed more succinctly as: 
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HF represents the transpose of the matrix [H__]. The 
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elements of ale may be evaluated by direct substitution of 
the expressions for the normal modes (Eqs. II.A.1 and II.A.4) 


into the definition (II.B.2). From Appendix B we write: 
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The scattering matrix for a single junction is defined in 
Appendix C to have matrix elements (which are themselves 


matrices) as follows: 


S S 
ey et fe Al2 4 


A2l A22 
z TyII sega et 
Sail 7 [vy +H H]~*Cyt-ty*7HI 
(II.B.3) 
. Tyllay-lpytyity - z 
Sa12 2ty' +H fl a 1 = AY = "S4448 
a -1,Tytt wf _ 
Sany = 2C1 + H(y7) 17" = HUI + Sail 
: ylig hgh yl ge ley eee ieee 
Byog = LI + acyl) ttyl ty-leacyty-dyty!? _ 7] = -H[S, 41H 


7) ere = 
iwi . Mp gh = 
is 7 hy ¥, ~ 


ar, ak «it sera Ap weseeaent ats eo sesces TC) KI 4 
| 7. vo snkiiens aie vette vd bersuteve- sé Ven ta ce re neo 

a teat haw bpK. 1 sepa: eibom iemron std es 19x 
Mette on a xibneqgh sort oe nodsiniret i 


= 
oo he | 
ay z 
Bat ale mi £03) tees! x* 7 + 
Ry: 


ntitcein: olgnts s as alban ee heaeiil - 


pate 
 & 


e_ i 
_ - 


a sshd! ainecele pb cyinet ever oF 5 tbe 


tevoltot 2s, (esos 


vt . 
Piz ao é 
" [ nil Lia. ] . {82 
| Ar. sth’ 
ae ae 
3 | 
‘iii oe rn + r= ie 
ea #f- a Phin fates + *yIs: = era® 
‘ ere gy (eR + tit = ..,285 


sar - 


* 


Note that each of the individual elements, H, Y, T, and I are 
themselves matrices. For a single fin-line strip there is 

a scattering matrix associated with the left boundary 
(minimum Z), [s,d, described above, a transmission matrix 
associated with the width of the strip, [T,J, and a 
scattering matrix associated with the right boundard, [S,J. 
The expression for an equivalent eee for the 
net effect of [s,J, [T.WJ, and [s,J, is developed in detail 


in Appendix D. The result is: 
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The equivalent scattering matrix for a fin-line with 
and arbitrary number of fins may be determined by repeated 
application of equations (II.B.4), with the appropriate 


transmission matrices used in the cascade. 
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III. COMPUTER AIDED DESIGN PROGRAM 


The CAD (computer aided design) program has two = eS 
1) design of a filter (description of physical dimensions) 
based upon user specified characteristics; and 2) analysis 
of a physically described fin-line filter. The user has 
the option of directly entering the analysis portion of the 
CAD if filter dimensions have been determined by some other 
means. In most cases, however, the user will specify the 
desired response of the fin-line filter (bandedge frequencies, 
passband ripple and upper skirt rejection) and the CAD will 
determine the filter dimensions. These dimensions are then 
automatically passed to the analysis portion of the CAD. 

The design portion of the CAD is built upon the 
theoretical work of Levy [Refs. 6, 7] and contains a modified 
form of a subroutine developed by Shih in [Ref. 2]. The 
design technique is fundamentally a conventional network 
synthesis. The insertion-loss characteristic of the 


distributed filter is approximated by the equation: 
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h is the passband ripple 


A_ is the guide wavelength 
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. is the guide wavelength at the filter center frequency . 
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AoA go? is the guide wavelength at filter bandedge, f 


So 


a is the first-kind Chebyshev polynomial of degree n 


nis the number of filter resonators. 


The insertion loss equation is evaluated for the input 
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Specifications to determine the number of resonators required. 


To properly determine the number of resonators for the general 


case, it was necessary to develop subroutines to solve 
Simultaneously the transcendental equations for the wave 
numbers of the eigenmodes. The functional dependency of t 
insertion loss equation, above, on guide wavelength is 
ambiguous for the Class II and Class III filter. Thus, it 
was alsc necessary to determine an "effective" guide 
wavelength which accounts for the fact that the wavelength 
changes as the wave propagates from one dielectric region 
into another. It is this effective guide wavelength which 
is used in Levy's formula. The required reflection 
coefficients of the n+l fins are then derived from the 
Chebyschev polynomial. Next, the scattering matrix for a 


Finite septum (Eqs. II.B.3) is evaluated by a bisection 


technique to determine the ntl fin dimensions which satisfy 
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these reflection coefficient requirements. Finally, the 
resonator dimensions are determined so that the electrical 
length of each resonator is (1/2) A so" 
With the dimensions of the filter thus determined, the 
CAD then performs a complete frequency response analysis 
based upon the mode-matching technique described in 
Beetions II.A and II.B. Tt should be noted that the equal- 
ripple Chebyshev polynomial expression very accurately 
matches the computed response when the filter is narrow band 
(less then 4% bandwidth). The agreement between the 
Chebyshev filter approximation and the mode-matching 
scattering analysis deteriorates as the filter bandwidth 
increases. A thorough discussion of this phenomenon is 
given by Shih in [Ref. 8] along with new formulae which 
account for the frequency dependencies of the step-impedances. 
However, excellent agreement between design criteria and 
predicted response has been obtained by employing an 
empirically derived offset to the lower bandedge frequency. 
That is, if the user specifies the filter bandedges to be 
Fl and F2, the design portion is actually computed for 


bandedge frequencies: 


a -= PiCl «= (-ssg=y—? Se and. F2. 


andar: 
vPro wa . 


The program listing given in Appendix F includes those 
subroutines which were developed as a consequence of the 
previously discussed theory. Subroutine TREQ develops 
eigenmode wave numbers for all three filter classes. 
Subroutine PROPAG determines the corresponding propagation 
numbers. Subroutine SMATRI determines the values of the 
H matrix coupling coefficients, and the elements of the 
scattering matrix. Subroutines CASCAD, DSUM, and TRANSM 
are used to determine the net effect of several cascaded fins 
and resonators, as discussed in Appendix E. 

Certain subroutines used in the design portion of the 
program are substantially the work of Shih LRef. 2), and are 
not included in the program listing. Also omitted are 


certain mathematical and algebraic routines required for 


the matrix manipulations. 


IV. PROGRAM VALIDATION 


A. OVERVIEW 

The first phase of program validation was conducted as 
critical subroutines became available. For example, the 
subroutine which determines wave mode numbers was tested 
for special circumstances for which the values could be 
verified by means of a hand calculator. A series of tests 
were then made to assure convergence to this limiting case. 
Professor Y. C. Shih kindly modified subroutines extracted 
from his work [Ref. 2] to output values at intermediate 
stages of his analysis routine. These numerical values 
were of interest because they were determined from closed- 
form expressions for the reflection coefficients, and 
employed a residue calculus technique. Thus, the corres- 
ponding values determined by the new program were evaluated 
from a completely different set of expressions and based 
upon a fundamentally different approach to the problem. 
Correspondence between equivalent values was considered 
good to within the accuracy of the numerical process. 

Predicted filter responses were compared to actual 
filter performance to validate further the computer aided 
design and analysis program. Dimensions for Class I and 
Class II filters described by Shih [Ref. 2] end Arndt, 
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the CAD and filter performance was predicted over an 
appropriate frequency range. The agreement between 
predicted response and published filter performance was 
excellent. Lastly, final validation of the CAD was accom- 
plished by means of locally designing, constructing and 


testing three Class I and two Class II fin-line filters. 


B. FILTER CONSTRUCTION AND TEST 

The locally constructed test filters were designed to 
be resonant within the 8 to 12 GHz band. This frequency 
band choice was selected to simplify both filter construction 
and test. 

The three Class I filters were manually cut with scissors 
and razor blades from a 2 mil(0.002 inch) .copper foil. The 
narrow band Class III filter was made by manually cutting 
away the 2 mil copper layer from 1/8" copper clad substrate 
(Rexolite 1422). The substrate for the wideband Class cre 
filter was a slab of styrofoam cut from a packing form. 
Conducting fins were cut from 2 mil copper tape and pressed 
onto the styrofoam slabe to form the filter structure. 

The test fixture used was a specially machined ten inch 
aluminum section of WR-90 compatible waveguide (0.9 by 0.4 
inch, internal), modeled after a design by Knorr LRef. ee 
This guide consisted to two longitudinal sections which 


were bolted together. The transverse dimension of the 
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Class I filter was made equal to the test fixture exterior 
It was a simple matter, therefore, to separate 


dimension. 
slightly the test guide halves, slip in the fin-line 
While this 
in the center area 


structure, and tighten the clamping bolts. 
were noted. 


filter was 


introduced a 2 mil dimension distortion 
Pie 
A small spot of 


of the test fixture, no adverse effects 
The fin-line structure of the Class 


pressed into the opened waveguide groove. 
rubber cement served to hold the dielectric substrate 


against the guide wall. 
Filter electrical performance was measured with the aid 


of an HP-8409 automated network analysis system. 


SULTS 
through 


Cx 
A plot of normalized throughput power versus frequency 
for seven filters is shown in Figures IV-1 through IV-7, 
Filters #1 


these 
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for Filters #1 through #7, respectively. 
fabricated by the authors. For 
plot of the predicted filte 


See | 
The circles contain 


#4 and #7 were 
CAD. 
Filters #5 and #6 are representa- 


the continuous 
response generated by the 


curve 

measurement data points. 
tive filters described by Arndt, et al [Ref. 3]. Again, 
the line is the CAD generated plot, but the data points 
LEtSErS . 


these curves were taken from figures contained in [Ref 


Table IV-1 contains -the physical descriptors of the 
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Filters #1 and #2 are Class III type structures; filters 
#3, #4 and #7 are Class I; filters #5 and #6 are Class II. 
In general, the shape of the filter response corresponds 
well to the simulated filter. Actual filter passband ripple 
was more than predicted, and varied from a best case of 
o) dB (#1) to a worst case of 3.0 dB C#2).., Best case 
insertion loss (Filters #1 and #7) was about 0.3 dB. 
Modifications to existing filters and multiple construction 
attempts at the same design suggests that the majority of 
the discrepancy in predicted and actual response was due to 
error in the hand-out filter dimensions. 

The good agreement between predicted and actual filter 
response is readily seen for filters #4, #6 and #7 (Figures 
IV-4+, IV-65, IV-7). In other instances the agreement is less 
apparent because of the expanded scale of the plot. For 
example, the discrepancy between predicted and actual 
response of filter #1 may be almost entirely accounted for 
by a shift in center frequency of about 150 MHz. This 
150 MHz shift represents about a 1.7% error in center 
frequency and could result from a construction error in 
resonator dimension of about 0.3 mm. 

Filter #7 is particularly interesting because it 
demonstrates the utility of the program. This filter was 
designed to satisfy a requirement of a physics researcher. 


Total time logged on the computer was about 15 minutes. 
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A copy of the output to the CRT is given in Appendix G&G. 

The CAD determined filter dimensions were scribed onto a 
sheet of copper foil and the filter was cut out by hand. 
Total construction time was about one hour. Figure IV-8 

is a copy of the test equipment output, with notations added. 
The minimum measured insertion loss was 0.28 dB at 8.604 GHz, 
Passband ripple was measured at 0.39 dB. The lower 3 dB 
point was measured at 8.360 GHZ and the higher 3 dB point 
was measured at 8.823 GHz. Thus, the CAD analysis values 

for the 3 dB frequency cutoffs of about 8.37 GHz and 

8.85 GHz compare to within 0.3% of the measured values. The 
researcher considered the more than 60 dB of out-of-band 


rejection provided by this filter to be quite adequate. 


ee. DISCUSSION 

Experiences related to the design, fabrication and 
test of the fin-line filters motivates us to catalogue a 
number of observations which may be useful to others. 

First, it should be stressed that the filters constructed 
were very crude by industry standards. Dimensional tolerances 
ere estimated at + 0.4 mm. Thus, the dimensional accuracy of 
the smallest fins is roughly bounded by + 40%. The reasonably 
good electrical performance of these filters was, therefore, 
all the more striking. Also, it should be noted that the 
approximate 150 MHZ shirt in center frequency for filter #1 
(narrow band Class III) can be the result of about a 0.2 mm 


error in the construction of the resonator dimension. 
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TABLE IV-1 


DESIGN DATA FOR CAD VALIDATION TEST FILTERS 

Filter # 1 2 3 - 5 6 7 
a (cm) 1.143 1.143 1.143 1,143 . 7895 . 188 1.143 
d (cm) 0.508 0.508 0.508 0.508 0127 -00635 0.508 
€,, (upper) 2, 550 1.030 1.000 1.000 26 12 1.000 
e, (lower) 1.000 1.000 2.000 2.000 2J22 22e 1.000 
Fl (GHZ) 8,625 9.950 8.450 11.30 15.00 66.0 8.400 
F2 (GHZ) Butte “£0595 8.800 11,75 15.60 67.0 8.800 
# Strips 3 6 4 4 4 4 7 
# Res 2 5 3 3 3 3 6 
Sl (cm) -0959 0.1485 0732 -4995 22079 - 0347 - 0906 
S2 (cm) -4928 0.7782 .4185 -1,4252 .7760 - 2038 .4794 
S3 (cm) 0959 1.0040 -4185 1.4252 .-7760 - 2038 5961 
S4 (cm) --- 1.0040 0732 4995 .2075 . 0347 -6134 
S5 (cm) --- 0.7782 --- --- ao “== .5961 
Sé (cm) --- 0.1485 --- --- --- --- .4794 
S7 (cm) --- --- --- --- --- --- - 0906 

(cm) 1.5044 1.1435 1.9991 - 9398 . 8400 -1895 2.0480 
R2 (cm) 1.5044 1.1142 2.0728 .9265 - 8400 sL900: 21227 
R3 (cm) ——— 1, LL04 <3. 5991 .9398 -8400 1895 2.1280 
R4 (cm) --- 1.1142 --- --- --- --- 2.1280 
R5 (cm) w-- 1.1435 --- --- --- --- 2eaaey 
R6 (cm) --- =e one nee ws aed 2.0480 
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Secondly, difficulties associated with the fabrication of 
the Class III filters suggests the need for an improved 
construction technique. Any new such technique needs to 
address the two major shortcomings of the method employed. 
These are: 

1) Difficulty in longitudinally aligning the two 
separate fin-line sections to maintain symmetry about the 
x=0 axis. 

2) Persistent problems associated with inadequate 
electrical contact between the filter's conducting fins and 
the waveguide's broad walls. 

Failure to correct either of these recurring tendencies 
resulted in severely reduced electrical performance. Thus, 
the original motivation for developing Class III filters 
may have been false. 

Thirdly, dielectric substrates should be tapered. For 
example, the narrow band Class III filter (filter #1) was 
originally constructed on rectangular dielectric slabs. 

The entire insertion loss curve of this filter exhibited 
about 12 dB of ripple. The insertion loss data points 
plotted are for the same filter, but with the addition 

of dielectric ramps on all four slab ends. Passband ripple 
was thus reduced to about 0.1 dB, and ripple along the 


rejection skirts was no longer noticeable. 
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correlation to passband ripple. That is, as the ripple is 
decreased, the fin-line strip dimension decreases. Strip 
width may be increased by either specifying a larger ripple, 
or by forcing additional resonators by requiring sharper 

rejection skirts. Wider strip widths are preferrable simply 
because of ee eee mechanical tolerances, and structural 


integrity. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The numerical analysis technique employed by the authors 
will always produce a result. The question is whether this 
result has any correlation to the physical world. It was 
found that this numerical analysis yields values essentially 
identical to those obtained by Shih [Ref. 2], who employed a 
completely different approach, and evaluated entirely dif- 
ferent expressions. This suggests that both approaches are 
valid. Furthermore, the good agreement between all 
predicted filter responses and filter performance strongly 
Suggests that the procedure applied was theoretically sound, 
and properly executed. 

The CAD program produced by this study is not an end in 
itself, but a design tool. It should be used to answer a 
number of design questions. It is recommended that an 
organized study be conducted to determine the following: 

1) Under what conditions is the center-line structure 
preferrable to a bilateral structure? 

2) For the bilateral fin-line, what determines the 
optimal transverse fin-line location? 

3) Are there optimal values for the relative permittivities 
in the Class II and Class Iil filter designs? 

4) How sensitive Se Filter reponse to dimensional 


tolerances? 
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5) Does the current filter design meet engineering 
requirements, or is the use of a design/response optimization 
routine justified? 

6) To what extent can the number of modes retained in 
the numerical analysis be reduced without adversely 
reducing accuracy? 

Finally, to determine more accurately the degree of 
correspondence between simulated and physical filter 
performance it will be necessary to construct higher quality 
devices. In particular, the fin-line structure should be 
produced by photo-etching techniques to minimize 


dimensional errors. 
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APPENDIX A 


REGION I EIGENMODES 


The solution to the general wave equation for an EM 
wave in a guide with disimilar dielectric zones, as shown 
in Fig. II-1 (Region I), is found in a manner similar to 
that for a simple waveguide. However, in this instance, 
the eigenvalues for the propagation numbers are not 
explicitly expressible. Instead, a pair of transcendental 
equations is found (A.1 and A.2) which must be solved 
Simultaneously to determine each eigenvalue. 

For incident TE Oo (m odd) modes, the Z component of the 
magnetic field (H_) in the wave guide will have a value of 
zero along the x=0 plane. It is convenient, therefore, to 
impose the condition of a "magnetic wall" along x=0, on 


which a boundary condition is that 


By referring to the geometry of the waveguide depicted 


in Fig. II-2, we can write the form of the E,, fields as: 
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At the magnetic wall, x=0, the boundary condition imposed is: 


"E: 
C-jwu Ae = — = — — oo 
Y x=0 
But since (s) = 0 for all fields, we can write: 
par) =0-= [KSA cosk. = k.. B sink: x1] : 
aX 7. lm m im im m a 7 
x=0 x=0 


Whence, A_=0. 
m 


The second boundary condition to be satisfied is that: 
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a that: 
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Therefore, the ES Fields in Region I is given by: 
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Or, employing a trigonometric identity, Ey may be expressed 
as: 
sin kom (a-x) 
6 = Cm __ = 2. Sin k__ (a-x) 
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Finally, Ey = Eu at the boundary of the two zones. That is: 
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Likewise, the Z component of the magnetic field must be 


continuous across the boundary. So that: 
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Substituting for Ba above, we find that: 
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The second equation which the propagation numbers must 
satisfy is the usual dispersion relation for an EM wave 


propagation through a dielectric discontinuity [Ref. 4]. 
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The amplitudes of the eigenmodes may be determined by 
applying the orthogonality property of the modes. That is, 
trigonometric identities may be used to rewrite the 


expressions for the eigenmodes as follows: 
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The orthogonality ofthese modes is exploited by integrating 


the product of mode pairs across the width of the guide. 
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APPENDIX B 


DERIVATION OF EXPRESSIONS FOR H-MATRIX ELEMENTS 


The elements of the H-Matrix are defined by: 
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APPENDIX C 


JUNCTION SCATTERING MATRIX 


Consider the waveguide discontinuity as a two-port junc- 
tion problem, as represented schematically in Fig. II-3, 


= + : 

where (a a a) and (b_ , b_) are the amplitudes of the 
n n n n 

incident and scattered modal fields in Region I and IT, 

respectively. 


The electric field is a single-valued function of Z. 


Hence, at the junction boundary, Z=0, we can write: 
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similariy; the relationship for the x-component of the 
magnetic field, Hs ge Be 
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Toutilize the orthonormality of the eigenmodes, we multiply 


Pal by - and integrate across the waveguide. 


We define the values a as: 
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Likewise, multiplying equation C.2 by ‘.. and integrating 


yields: 


The wave impedances of the normal modes, as may be repre- 
sented as diagonal matrices, the values _ ee as anm by n 
matrix, and the amplitudes of the normal modes as vectors. 


Equations C.3 and C.4 may then be expressed as follows 


(subscripts have been omitted for simplicity). 
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But from standard scattering matrix notation with the 


values denoted as in Fig. II-3, we can also write: 
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APPENDIX D 


SCATTERING MATRIX FOR A FINITE LENGTH SEPTUM 


A generalized scattering matrix technique in a manner 
similar to that of Shih in [Ref. 3] is applied to obtain the 
two port scattering matrix for the finite septum. 

Consider a ce -type wave incident from Region I at the 
junction A. The field at this junction is comprised of the 
vector sum of all the eigenmodes in each region. Some of 
the energy will be reflected at the boundary, propagating 
in ae modal fields in the minus Z direction. The remainder 
of the energy will be transmitted into Region II in oe modal 
Fields in the positive Z direction. After traveling a 
distance D, the transmitted modes encounter the second 
junction, B, where again some energy is reflected, and some 
transmitted. Multiple reflections of the modes between 
junctions A and B is inherently accounted for in the matrix 
manipulations which yield the scattering matrix for the 
finite septum. 

In C a single junction was analyzed as a two-port scat- 
tering matrix, [sJ. Each of the four elements of this 
scattering matrix, [S Toms bese ES J, ana [s 
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a matrix of infinite dimensions, corresponding to the 
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infinite set of eigenmodes in both regions (equation C.8). 
The junction B is identical to A, except for orientation. 


It may be shown, therefore, that the scattering matrix for 
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The effects of the wave propagating (for propagating modes ) 
or attenuating (for evanescent modes) along the strip 
length are accounted for by a transmission matrix, [TWJ,; 


defined as: 


where [I] is the identity matrix, [9] is the null matrix, 
and (T] is a diagonal matrix of infinite size. The 


diagonal elements of [T] are: 
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The combined effect of the single junction A followed by the 


transmission along the strip for a distance D is given by 


the scattering matrix ist, J. 
’ = iE a 
[Ss ~ (TJ oe (TJ 
Finally, the net effect of the finite septum is determined 
as follows. 
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With the equalities given in equations D.5 and D.6, we 


substitute equation D.2 into D.3 to yield: 
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The matrices which are the elements of the equivalent 
scattering matrix may be written down by simply comparing 


equations D.9 and D.10 with equations D.7 and D.8. 
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APPENDIX E 


CASCADING SEVERAL FINITE LENGTH SEPTA 


The fin-line filter is simply a series of finite length 
septa separated by transmission line sections. The net 
effect of the entire filter structure may be represented in 
a single "total" scattering matrix which accounts for the 
combined effects of each septa/transmission line section. 
This total scattering matrix is obtained by combining each 
septum equivalent matrix in a manner precisely analogous 
to that used. in D to develop the equivalent matrix fora 
finite septum. That is, in D, expressions for the net 
effect of two discontinuities separated by a transmission 
line section were obtained (equation D.11). It is now 
necessary to express the net effect of two septa separated 
by a transmission line section. Consequently, the equiva- 
lent scattering matrix of two separated septa may be found 
by use of the equation D.11, with the appropriate 
scattering matrices for the septa substituted for the 
junction scattering matrix. More than two septa are 
cascaded by repeated application of D.11. Due care must be 
exercised to assure that the propagation numbers used in 
the transmission matrix correspond tc the appropriate wave 


number and region of waveguide propagation. 
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